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Abstract 

We attempt to see how closely we can formally obtain the planetary and 
light path equations of General Relativity by employing certain operations on the 
familiar Newtonian equation. This article is intended neither as an alternative 
to nor as a tool for grasping Einstein's General Relativity. Though the exercise is 
understandable by readers at large, it is especially recommended to the teachers 
of Relativity for an appreciative understanding of its peculiarity as well as its 
pedagogical value in the teaching of differential equations. 



Everyone knows Newton's theory of gravity and some know Einstein's the- 
ory of General Relativity (GR) . Undoubtedly, GR is one of the most beautiful 
self-consistent modern creations in the realm of theoretical physics. It has won- 
derfully tested against various astronomical observations to date including those 
in the Solar system. However, at a popular level, a naive question is often asked 
as to whether the GR effects could have been interpreted using a more mundane 
theory than the abstract theory of GR in which gravity - which is as real a force 
as any other - has been "geometrized" . For instance, some ask the question: 
What is the difference between the bending of light rays in GR with that occur- 
ing in a refractive optical medium? The answer lies in the well known fact that 



1 



the propagation of light rays in a gravity field a la GR can be exactly rephrased 
as propagation in an equivalent optical refractive medium with appropriate con- 
stitutive equations [1]. The refractive index can be employed in a new set of 
optical-mechanical equations so that a single equation covers motions of both 
massive and massless particles in a spherically symmetric field [2-4]. An ap- 
proach of this kind provides a useful and interesting window to look at familiar 
observed GR results but, by no means, implies a replacement of GR. 

The whole point of the above paragraph is that one inevitably needs to know 
the metric solutions of GR in advance. Only after knowing them, one can derive 
appropriate refractive indices and the method of optical-mechanical analogy in 
terms of these indices then exactly reproduces the GR geodesic equations. That 
is to say, we might employ different working methods but the physical content 
remains essentially that of GR. (There have been attempts to set aside GR 
altogether and propose alternative physics by introducing a variable test mass 
[5] , or even assuming variable speeds of light in flat space [6] . These ideas have 
their own values and we are not going to discuss them here.) 

The object of the present article is somewhat different: We are not going 
to suggest any working method of the kind described above, but present an 
interesting calculation. (However, it must not be weighed against the grand 
edifice of GR). Using PPN-like approximations on the Newtonian theory, we 
shall formally obtain planetary and light path equations. They resemble the 
path equations of GR only fortuitously and this is the amusing part. Apart from 
this, the contents could be instructive in exemplifying the role of numerically 
smaller terms in the differential equations. 

To begin with, one recalls an earlier discussion of M011er [7] that has shown 
that the bending of light rays is due partly to the geometrical curvature of space 
and partly to the variation of light speed in a Newtonian potential. In fact, the 
ratio is exactly 50:50. The GR null trajectory equations can be integrated, once 
assuming a Euclidean space with a variable light speed and again a curved space 
with a constant light speed. This analysis and arguments clearly elucidate the 
complementary roles of curved space and Newtonian theory in the best possible 
manner. This complementarity motivates us to examine how far, if at all, we 
are able to introduce curvature effects in the path equations of the Newtonian 
theory. That is: We try to obtain, from the familiar Newtonian theory itself, the 
form of the known GR path equations of motion without geometrizing gravity. 
(It is known that the gravitational rcdshift is a prediction of GR, but it is also 
known that it can be predicted from the Equivalence Principle without using 
GR equations [8]. Hence we shall not address this result here.) 

Let us start from the usual Kepler problem of a massive test particle moving 
around a spherical gravitating mass M under the Newtonian inverse square 
law. Let T and V denote the kinetic and potential energies respectively. Then 
T + V = constant = (say) implies in relativistic units 

1 r' 2 , 2 ■ 2i 2-1-^0 / -i \ 

-\r +r if }-mc r = — (1) 

where m = GMc^ 2 and a dot denotes differentiation with respect to Newtonian 
time t, Co is the speed of light in vacuum. The central nature of the force implies 
constancy of the angular momentum (the Lagrangian is independent of <p) such 
that 

r 2 cp = h . (2) 
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With u = p, we can rewrite Eq.(l) as 

- 2mucl = E Q (3) 

where the constant Eq has the dimension of Cq. For bound material orbits Eq 
< 0. Customarily, by differentiating again with respect to <p, one finds a second 
order differential equation that yields a Keplcrian ellipse given by 

u = -(1 + ecosip) (4) 
P 

where e is the eccentricity, p = is the semi-latus rectum. 

Let us redefine the radial variable u — > u' through the equations 

u' = (5) 



hi 



duY 



*(«) = (1 + — 



r 



(6) 



(7) 



(Aside: These transformations are not unfamiliar to those conversant with GR.) 
After some straightforward algebra, we get 



where 



dv! = &(u)tt(u)du 
, w . / mu\- 1 / mu\ 



fi(u') = (1 - 2mu') 



'\2 



1 + (1 - 2mu')5 



(8) 

(9) 
(10) 

(11) 



Note that $(u) of Eq.(6) is numerically the same as $(m') of Eq.(ll). The same 
applies between Q(u) of Eq.(9) and Q(u') of Eq.(10). The following expansions 
can also be directly verified: 

2mu = 2mu' + 2m 2 u' 2 + 5m 3 u' 3 + ... = 2mu' + 0(m 2 u' 2 ). (12) 

This implies that, to first order, r ~ r'. Also, 

1> 2 (u')n 2 (u') = 1 - Amu' + 0(m 2 u' 2 ). (13) 

Let us now express Eq.(3) in terms of the new variable vl . Multiplying both 
sides of Eq.(3) by $ 2 f2 2 and using Eqs.(5)-(13), we get 



il u + 



dv/_ 

d(p 



= c 2 [E Cq 2 + 2mu' + 0(m 2 u' 2 )] $ 2 ft 2 . (14) 
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Simplifying further using Eqs.(lO) and (13), we have 



hi 



u' 2 +{^) -2mu' 3 
dip 



= c 2 [E q Cq 2 + 2mu'{\ - 2E Co 2 ) + C>(mV 2 )] 



(15) 

Apply this equation to a practical situation, the Solar system. At the site of 
Mercury, the planet nearest to the Sun, mu ~ mu' ~ 2.5 x 1CT 8 . Let us ignore 
the terms 0(m 2 u' 2 ) in comparison to the mu' term. Then Eq.(15) reduces to 



h 2 



, 9 / du' 
u' 2 + 



V d<p 



2mW 



E a + 2mu'c 2 Q (l~2E c^ 2 ). 



(16) 



Differentiating with respect to <p, we get 



u' + ^ = ^- + 3mw 

dip^ p' 
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(17) 



where 



1 
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ho 



h' 2 



(l-2E oCo - 2 )i 



(18) 



is a rescaled constant. 

The final Eq.(17) seems suggestive with the usual perturbation term 3mu' 2 
appearing: It is exactly of the same form as the GR path equation! One notes 
that the constant h' involves the test particle energy E$ similar to what one finds 
in the GR treatment. To see this, compare with Eq.(17) the corresponding GR 
equation given by (Take henceforth G = 1): 



u + 



d 2 u 1 
d^ 2 



+ 3mw 



(19) 



in which p is given by p 



Mmicf. 



where m is the test particle rest mass, 



J = — is the constant angular momentum rescaled by the energy at infinity 
Ua = i m °^° and the constant U3 = r 2 ^-, A being the affine parameter [91. 

V 1 " r oo/ C 

As usual, considering low velocity, we can take J7o = m c§ and identifying the 
asymptotic value of J as ho, we have 



h 2 

y M 



(20) 



With this value of p, the GR perturbation term imu 2 then gives the well known 
perihelion advance of the Keplerian ellipse. 

In our case, the parallel of p from Eq.(17) is: 



p := 



meg 



hi 



M(l - 2E Q c ■ 



(21) 



Its asymptotic value can be computed using Eq.(l). For near circular orbits, 
the kinetic and potential energies are roughly of the same order of magnitude 
such that the velocity is v 2 <~ = muc 2 ,. Then, from Eq.(l), and noting that 
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u~u' asymptotically, we can write E — amu'c^ where — 1 < a < 1. Then the 
denominator becomes M(l — 2amu'). The term 2amu' ~ 1CP 8 can be easiliy 
ignored compared to unity and we are left with 

'=4 < 22 » 

just as in Eq.(20). So we can replace p' in Eq.(17) by its asymptotic value p 
given either by Eq.(20) or (22). 

For the motion of light, the situation is different: the dimcnsionless quantity 
EqCq 2 must be fixed to the value ^ so that p' — > oo. Recall that only a nonzero 
value for light (E$ ^ 0) in Newtonian theory is consistent with the zero value in 
GR [10]. (The zero rest mass of photons is a Special Relativistic or GR concept 
but is not a Newtonian concept). Consequently, we have the equation of the 
light ray trajectory exactly as in GR: 

v! + ^4r = 3mu' 2 (23) 
dp 2 - 

Thus Eqs.(17) and (23), respectively, seem to provide the same GR results 
as far as the weak field tests for the perihelion advance and the bending of light 
are concerned. To examine the situation more closely, recall what steps were 
involved. The first step is the radial rescaling u — ► u' which has no physical 
import. The second step is that, in arriving at Eq.(16), we have ignored terms 
like 0(m 2 u' 2 ) on numerical grounds. Note that it is only Eq.(15) per se, and 
not Eq.(17), that inverts exactly to the original Eq.(3) in the (u, ip) coordinates 
describing the inverse square law. As we see, Eq.(17) produces an additional 
3mu' 2 term! Strictly speaking, Eq.(17) is approximate to the extent we ignored 
the smaller terms compared to unity (of the order of 1CT 16 and less!) in arriving 
at it. Treating this Eq.(17) as an exact equation means that we are retaining the 
cubic additional term as the only perturbation while disregarding the remaining 
smaller perturbations. This is the only nontrivial step we have adopted in the 
above computation. 

If we had retained the smaller terms in Eq.(15), then it could tell the original 
situation: the exact Newtonian orbits. It is our nontrivial, but numerically 
justified, omission of the smaller terms that has brought forth equations similar 
to those in GR. Thus the exact solution of Eq.(15) is still a Keplerian ellipse 
but its expression does not look as familiar as in Eq.(4). Instead, in the primed 
coordinates, it looks like 

v! = u$(u) = — — (1 + ecos<p). (24) 
P 

where u is given by Eq.(4). Expressions might differ in looks depending on the 
choice of coordinates, but the orbital shapes do not change. 

One might think that though Eq.(17) looks different from Eq.(15), it still 
represents a Keplerian ellipse in the (u',ip) coordinates. This is not the case 
since Eq.(17) is now nonlinear. We can find its solution by standard procedures 
starting with the zeroth order solution u' = i(l + e cos ip) which is the solution 

of v! + = |. Eq.(17) then gives the observed perihelion advance as ^p-. 
[Note that if one starts with the same u' in Eq.(15) or its second derivative 
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form, one would eventually end up with Eq.(24) as the final solution]. Likewise, 
the exact equation for a straight line is 

u' = -4*( u ) cos< P C 25 ) 
it 

where R is the distance from the origin. To zcroth order, u' = is a 

solution of v! + = 0. By usual methods again with Eq.(23), one finds a 
total observed bending of light rays Aip ~ ^jf. 

The procedure leading to Eq.(17) has some similarity with that in GR. In the 
curved spacetime of GR, one needs to consider coordinate independent proper 
length I instead of the radial coordinate r. Thus, in the Schwarzschild metric, I 
is given by 



I = 



dr 



^/r(—2m + r) + 2m^/2m — r arctan yjrj (2m — r) 

= (26) 



'-) 



In terms of (l 7 ip) coordinates, the GR Eq.(19) can not maintain its form or 
assume another exact closed form due to the fact that r can not be expressed in 
terms of / in a closed form. However, in the weak field region, r ~ I, and we can 
maintain the form of Eq.(19) as it is, while ignoring higher order terms in I. In 
the present calculation, the background is Euclidean and so we can express I, 
using Eq.(8), as 1= J dr— J <J>(r')f2 _1 (r')dr' . In our calculation, we have ignored 
higher order terms in u'in the weak field region so that r ~ r'and we ended up 
with Eq.(17). 

Can we physically interpret our nontrivial step as a modification of the New- 
tonian force law? In this context, it is to be noted that, historically, Newton 
himself attempted to modify his force law to explain some phenomenon (for 
details, see Ref. [5]). One might also recall other efforts, for instance, Sommer- 
feld's calculation [11] for the precession of an electron in a Coulomb potential 
due to a proton (Z = 1): 

d ( m lt \ __Ze> ? 



dt \VT~i 

where r is a unit vector in the radial direction and e is the electronic charge. 
However, it produces only (l/6)th of the observed perihelion advance of planets 
if the Coulomb potential on the right is replaced by the Newtonian potential. 
One could try the above special relativistic equation with another kind of force 
law on the right [12] 

d ( m Q ~v \ Mm _ 



dt \^l-v 2 J r 2 (l-v 2 )i 

r 2 + r 2 ip 2 does produce the observed perihelion advance, but the 
difficulty is that its first integral does not produce the conserved relativistic 
energy. This is understandable because the potential is velocity dependent. 



6 



Coming back to our calculation, one might say that Eq.(17) [which is the same 
as Eq.(19)] corresponds to a potential V(r) = —— — ^ but then the last term 
leads to a dimensional mismatch (sec ref.[5]). Because of this, our procedure 
can not be interpreted as a modification of the Newtonian force law. Also, there 
was absolutely no use of the concept of geometric curvature in the calculation; 
it was completely Euclidean. 

Thus, we conclude that the similarity between Eqs.(17) and (19) is only a 
fortuitous though amusing coincidence; it is just a mirage resulting from the 
choice of coordinates. There is absolutely no reason to prefer (u',tp) coordinates 
over others and in this case, the formal coincidence will be lost. Nonetheless, 
the procedure illustrates something of pedagogical importance in the treatment 
of differential equations: One should be watchful with smaller terms! Their 
removal can nonlinearize a given linear equation [like going from Eq.(15) to 
(17)] and conversely, their restoration can linearize a known nonlinear equation 
[like returning from Eq.(17) to (15)]. 

It is a pleasure to thank Guzcl Kutdusova and Arunava Bhadra for useful 
discussions. 
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